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Abstract 

An  analytical  study  is  presented  regarding  the  theoretical 
behavior  of  an  airfoil  in  a  pitching  airflow.  The  theoretical  develop¬ 
ment  includes  the  momentum- integral  equation  for  boundary  layers  in 
unsteady  flow  and  the  von  Karman-Pohlhausen  integral  method  for 
unsteady  flow.  A  computer  program  was  written  to  model  this  and  was 
applied  to  a  symmetric  Joukowski  airfoil.  The  results  of  this  study 
show  a  linear  relationship  between  the  increase  in  the  stall  angle  of 
attack  and  the  non-dimensional  pitch  rate,  ^Cu/I^. 
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I  Introduction 


0 


Discussion 

Since  the  earliest  days,  aerodynamicists  have  been  interested  in 
the  performance  of  a  wing  under  dynamic  angle  of  attack  changes.  The 
earliest  work  seems  to  have  been  most  concerned  with  these  effects  as 
they  pertained  to  an  aircraft  experiencing  rapid  changes  in  angle  of 
attack  due  to  the  encountering  of  a  gust  (Ref.  1).  More  recent  work 
seems  to  be  concerned  with  these  effects  as  they  pertain  to  the  motion 
of  a  wing  in  otherwise  still  air,  and,  in  particular,  this  effect  as 
applied  to  oscillatory  motion  as  might  be  experienced  with  helicopter 
blades  (Ref.  10). 

Both  of  these  dynamic  conditions  lead  to  an  enhanced  maximum 
lift  effect;  this  effect,  however,  can  be  separated  into  two  regimes. 

In  the  first,  the  coefficient  of  lift-vs-angle  of  attack  curve 
(C£-vs-a)  is  extended  up  and  to  the  right  so  that  the  maximum  lift 
coefficient  (Cimax)  is  increased  by  some  amount  ACJlmax.  In  both  kinds 
of  dynamic  angle  of  attack  change  cases,  the  Cimax  seems  to  be 
directly  relatable  to  the  non-dimensional  angular  rate  %Ca/Uoo,  where  C 
is  the  airfoil  chord,  d  is  the  pitch  rate  of  the  gust  or  the  wing,  and 
Uoo  is  the  freestream  velocity.  In  the  second  regime,  a  large  ACfmax 
appears  as  a  spike  superimposed  on  the  C£-vs-a  curve  followed  by  a 
catastrophic  decline  in  C i.  This  later  effect  has  been  attributed  to 
the  appearance  of  strong  leading  edge  vortices  (Ref.  10),  and  is  beyond 
the  scope  of  this  effort. 

This  study  is  concerned  with  only  the  first  regime.  While  the 
general  trend  of  the  effect  in  this  regime  is  similar  under  both  motions 
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there  appears  to  be  a  measurable  difference  in  the  extent  of  the  effect 
even  if  the  restriction  to  constant  angle  of  attack  change  rate  (a)  is 
fixed.  In  the  case  of  a  gust,  constant  a  means  that  che  wing  is 
stationary,  and  the  flow  field  is  made  to  change  direction  at  a  constant 
rate.  While  in  the  case  of  the  moving  wing,  the  flow  is  fixed,  and 
the  wing  pitches  into  it  at  a  constant  rate  from  some  pivot  location 
commonly  taken  as  the  half  chord  (Ref.  9). 

Theoretically  speaking,  these  two  conditions  differ  in  frame  of 
reference,  the  latter  being  more  difficult  to  analyze.  This  study 
is  limited  to  the  theoretical  examination  of  the  first  condition;  that 
of  a  fixed  wing  encountering  a  gust. 

Since  the  phenomena  of  stall  is  due  to  separation  of  the  flow  from 
the  upper  surface  of  the  wing,  it  is  clear  that  a  study  of  dynamic 
stall  must  include  and,  in  fact,  hinge  upon  a  study  of  the  boundary 
layer.  The  study  is  limited  to  laminar  flow,  but  tie  techniques  should 
be  applicable  to  turbulent  boundary  layer  flow. 

Problem  Statement 

To  set  up  the  problem,  the  momentum- integral  equation  for  boun¬ 
dary  layer  in  steady  flow  is  extended  to  allow  for  unsteady  flow. 
Following  the  theoretical  development,  a  computer  program  will  be 
developed  to  allow  one  to  apply  this  theory  to  certain  specialized 
problems;  namely,  symmetric  airfoils.  The  application  in  this  study 
was  to  a  Joukowski  J011,  an  11%  thick  symmetric  airfoil,  in  order  to 
model  the  flow  about  the  leading  edge  of  a  NACA  0012  airfoil,  popular 
in  helicopter  blade  use.  Finally,  the  theoretical  results  are  compared 
against  experimental  data  available  in  the  literature  (Ref.  1  and  9) 


Assumptions 

In  order  to  perform  the  analysis,  a  simplifying  assumption  that 
the  flow  is  incompressible  is  made.  In  addition,  since  the  point 
of  flow  separation  on  the  airfoil  and  the  relationship  to  be  loss 
of  lift  (aerodynamic  stall)  is  somewhat  arbitrary,  in  this  study, 
aerodynamic  stall  is  defined  as  when  flow  separation  (as  indicated  by 


the  shape  parameter)  occurs  at  the  quarter  chord  point  on  the  airfoil. 


II  Background 


In  canvassing  past  efforts  on  the  dynamic  stall  phenomenon,  it 
becomes  apparent  that  the  majority  of  the  work  concerns  the  oscillatory 
motion  on  an  airfoil,  which  is  most  closely  associated  with  helicopter 
blades.  Other  works  included  research  on  airfoils  pitching  into 
airflows.  Only  a  few  papers  covered  work  on  the  gust  problem,  i.e. 
airflow  pitching  around  an  airfoil.  The  genisis  of  this  work  was  an 
experiment  by  Max  Kramer  (Ref.  1). 

The  Kramer  Experiment 

In  1932,  Max  Kramer  published  a  paper  regarding  his  experiment 
on  the  dynamic  stall  of  airfoils  in  pitching  airflows.  He  got  the  idea 
after  listening  to  several  pilots  commenting  about  the  behavior  of 
wings  in  unsteady  airflows.  To  investigate  these  observations,  he 
set  up  an  experiment.  In  it  he  fixed  an  airfoil  between  a  diffuser 
and  a  set  of  movable  guide  vanes  (see  Fig  2) .  The  air  flowed  from 
a  source,  through  the  vanes,  over  the  airfoil  and  was  collected  in 
the  diffuser.  The  insertion  of  the  vanes  allowed  Kramer  to  alter  the 
direction  of  the  airflow  while  taking  data.  This  allowed  him  to 
simulate  a  gust  encountered  by  an  airfoil  in  otherwise  steady  flow. 
Kramer  conducted  his  experiment  using  three  different  airfoils — two 
different  profiles  of  the  same  chord  and  one  airfoil  of  the  same 
profile  as  one  of  the  others  but  having  a  different  chord.  Kramer  also 
ran  the  air  at  two  different  dynamic  pressures,  enabling  him  to  obtain 
the  pertinent  characteristics  involved  in  the  problem,  i.e.  effect  of 
velocity,  area,  camber,  etc.  To  correlate  the  data,  Kramer  used  a 
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non-dimensional  pitch  rate  parameter: 


where  C  is  the  chord  of  the  airfoil,  I^q  the  freestream  velocity  and  a  ' 
is  the  pitch  rate  of  the  airflow  in  radians  per  second. 

Based  on  his  experimental  results,  Kramer  developed  an  empirical 
relationship  between  the  increase  in  C£max  in  the  gust  over  C£,max  in 
steady  flow  and  pitch  rate 


C  da 


Cfcmaxdyn  55  c^maxst  +  0.36  —  ^ 


(Ref.  1) 


The  results  of  Kramer's  experiment  have  been  reproduced  in  Figure  3. 


The  Deekens  and  Kuebler  Experiment 

Another  experiment  of  interest  is  one  conducted  by  two  students 
at  the  Air  Force  Academy.  In  it,  an  airfoil  (NACA  0015)  was  emplaced 
in  a  low  speed  smoke  wind  tunnel.  Once  steady  flow  had  been  established, 
the  airfoil  was  pitched  at  a  constant  rate  through  a  servo  motor  con¬ 
trol  (see  Fig.  A).  As  in  the  Kramer  experiment,  a  direct  correlation 
was  found  between  the  non-dimensional  pitch  rate  and  the  delay  in 
stall  angle  of  attack.  As  stated  in  the  introduction,  this  is  a 
different  problem  than  that  of  the  rotating  airflow.  The  results  as 
shown  in  Fig.  5  can  be  expressed  as: 

Aostall  *  2.1  (^Co/Uqjj), 

where  C  is  the  chord,  a  the  angular  pitch  rate  and  the  freestream 

velocity.  Applying  classical  airfoil  theory,  this  translates  to: 

ACimax  =6.57  Co/Uqq,  (Ref.  9:2-15) 
a  much  higher  slope  than  that  obtained  by  Kramer. 
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This  is  a  clear  indication  that  there  are  different  physics 
governing  the  behavior  of  the  airflow  in  each  problem,  and  both  should  be 
pursued.  This  study  will  investigate  the  behavior  of  the  airflow  and 
boundary  layer  of  the  first  case. 


Ill  Theory  Development 


Introduction 

In  deriving  the  theory  for  the  gust  problem,  use  was  made  of  the 
procedures  used  to  derive  the  momentum- integral  equation  for  the  boun¬ 
dary  layer  in  steady  flow.  That,  derivation  is  reproduced  in 
Appendix  B. 

The  momentum- integral  equation  for  the  boundary  layer  is  another 
form  of  the  boundary  layer  equations  that  result  from  the  simplifications 
of  the  Navier-Stokes  equations  that  govern  all  flows.  Thus,  a  procedure 
similar  to  that  used  to  derive  the  boundary  layer  equations  in  general 
will  be  used  here;  that  is,  to  develop  the  continuity  equation  and 
the  momentum  equation  for  a  fluid  element  in  the  boundary  layer. 

The  Continuity  Equation 

In  words,  the  continuity  principle  states  that  the  net  rate  of 
mass  flow  out  of  a  control  volume  is  equal  to  the  time  rate  of  loss 
of  mass  within  that  control  volume  (Ref.  2:  10).  Referring  to  Fig.  7, 
this  gives: 

h  h  h  h 

/  pudy  +  _3_  (Zpudy)dx  -  /pudy  +  motop  =  -  3  /  pdydx 
o  3x  o  o  3t  o 

Reducing  this  equation  gives: 

h  h 

motop  *  -  _3_  /pdydx  -  _3_  (/pudy)dx  (1) 

i  3t  o  3xo 


The  Momentum  Equation 


The  momentum  principle  states  that  the  sum  of  the  forces  acting 
on  the  fluid  in  a  control  volume  is  equal  to  the  net  rate  of  transport 
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of  momentum  out  of  the  control  volume  plus  the  time  rate  of  change  of 
momentum  In  the  control  volume  (Ref.  2:11).  Referring  to  Fig.  8,  the 
x-component  (parallel  to  the  wall)  of  the  momentum  equation  may  be 
found  by  finding  the  sum  of  the  forces  in  the  x-direction: 

h  h  h 

£F  =  -xwdx  +  / pdy  -  /pdy  -  )  (/pdy)dx, 
o  o  3xo 

and  equating  it  to  the  momentum  terms 

h  h  h  h 

*  /puudy  +  _3_(/puddy)dx  -  /puudy  +  m^  •  Ue  +  _3_  /pudydx 
o  3x  o  o  3t  o 

such  that 


h  #  h  h 

3_  /(puu)dydx  +  m^  *  Ue  +  3_  /pudydx  =  -  xwdx  -  3_  / pdydx  (2) 

3xo  p  3t  o  3xo 

Substituting  for  m^  from  equation  (1),  equation  (2)  becomes: 

h  h  h  h 

3  /  (puu)dydx  -  Ue  *  3_  /pdydx  -  Ue  •  3_  /pudydx  +  3_  /pudydx 

3xo  3t  o  3xo  3t  o 


h 

*  -  xwdx  -  3_  /pdydx 
3x  o 

Thus,  for  this  two  dimensional  problem,  the  pressure  gradient 
in  the  x-direction  is  the  only  pressure  gradient,  and  correspondingly, 
the  partial  derivatives  of  pressure  are  taken  as  total  (not  substantial) 
derivatives.  Also,  since  x  and  y  are  independent,  the  partial  deriva¬ 
tive  with  respect  to  x  can  be  brought  across  the  integration  symbol, 
since  the  integration  is  in  the  y  direction.  The  equation  becomes 

h  h  h  h 

/  3  (puu)dydx  -  Ue  *  /pdydx  -  Ue  .  /  3  (pu)dydx  +  3_  /pudydx 
o  3x  3t  o  o  3x  3t  o 

h 

=  -  xwdx  -  /  d£  dydx 


r:-  ■ 

i-  ■, 


h-- 


ffl 


h 


n  o 


N 


If  we  further  restrict  the  flow  to  incompressible,  the  density  can  be 
divided  through,  as  can  the  common  dx  term,  so  that 


^  8  (uu)dy  3  8  udy  8  ^ 

/  —  -  Ue  .  ^  /d,  -  Ue  .  /  jj  +  / udy 


8x 


(4) 


p  dx 
o 


To  this  point,  the  only  difference  between  this  equation  for 
unsteady  flow,  and  the  equation  for  steady  flow  is  the  addition  of  the 
time  dependent  terms.  However,  the  next  step  is  one  of  the  key  develop¬ 
ments  in  the  momentum-integral  equation  for  the  boundary  layer  in 
unsteady  flow.  It  was  at  this  point  in  the  steady  flow  equation  that 
Euler's  equation  was  introduced  to  substitute  for  the  pressure  gradient, 
and  the  same  will  be  done  here.  As  shown  in  Appendix  A: 


_  I  .  Ue  +  3Ue 


(5)  and  (A13) 


p  dx  ~~  9x  8t 
The  additional  9Ue/8t  term  introduced  here  is  significant  and 
will  be  shown  to  be  a  contributing  factor  in  the  behavior  of  the  flow 
field  in  gusts. 

It  is  important  to  note  that  the  potential  flow  field  in  incom¬ 
pressible  flow  can  and  does  respond  instantaneously  to  changes,  such 
as  an  angle  of  attack  change.  This  is,  if  an  airfoil  is  inserted 
in  a  stream  at  some  angle  of  attack  aj.,  and  the  airfoil  is  allowed  to 
pitch  up  to  an  angle  of  attack  a2»  the  flow  field  will  be  the  same  as 
if  the  airfoil  was  placed  in  a  freestream  at  an  angle  of  attack  02  with 
no  pitching.  This  is  because,  as  shown  in  Appendix  A,  the  velocity 
in  potential  flow  is  a  function  of  a  only  and  is  not  a  function  of 
time.  However,  the  pressure  gradient  does  not  respond  in  this  fashion. 
The  flow  field  resists  changes  due  to  the  inertia  of  the  flow,  in 
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addition  to  the  time  dependent  term  introduced  here  which  causes  a 
different  kind  of  behavior  for  the  pressure  gradient. 

Continuing  on  with  the  derivation,  now  substitute  equation  (5) 
into  equation  (4)  to  obtain: 

r  9  (uu)dy  9  ^3  udy  a  ^ 

’  *  -  Ue  n  'dy  -  Ue  '  H  +  u  /udy 

o  o  o  o 

. .  I-  +  J0e  Jy  + )  m  dy 

0  0  X  0  L 

o  o 

which,  after  rearranging,  yields 

^  9  (uu)dy  ...  9  ...  ^9  udy  9  ^  , 

-  ’  3l  +  Ue  /dy  +  Ue  /  -  u  /udy 

o  o  o  o 

J.  M,  ,  *  SUe  .  T™.  (6) 

+  /Ue  dy  +  /  -jj  dy  - 

o  o 

Looking  only  at  the  time  dependent  terms 

Ue  h  'dy  +  '  ^ dy  -  It  'udy’ 

o  o  o 

note  that  the  first  term,  h,  is  an  arbitrarily  selected  constant.  When 
the  integration  is  carried  out  and  evaluated  at  the  limits  and  then  the 
partial  derivative  with  respect  to  time  is  taken,  the  term  is  zero. 

In  the  second  term,  since  the  y-coordinate  is  independent  of  time,  the 
partial  with  respect  to  time  can  be  taken  across  the  integration  sign 


so  that 


/  dy  -  /udy  =  |^  /(Ue  -  u)dy. 
o  o  o 


Now  examine  the  spatial  terms  of  the  left  hand  side  of  equation  (6): 


*!  8  (uu)dy  .  „  *!  8u  .  ,  ,,,  8Ue  . 

-  '  *  y  +  Ue  /  dy  /Ue  dy. 

o  o  o 


10 


« 

—  1 


t  ■ 

I.  ' 
l-' 


t 


The  second  term  of  equation  (8)  looks  like  part  of 

o 

so  that  expanding  this  term  out  yields 

h  a  h  »»»  h  - 

/  (uUe)dy  =  / u  “  dy  +  /Ue  —  dy. 
o  o  o 


(9) 


Note  that  the  second  term  on  the  right  hand  side  of  equation  (9)  can 
be  written 


/Ue  dy  -  /  dy  -  /u  dy. 


3 (uUe) 
3x 


3Ue 

3x 


(10) 


(ID 


Substituting  (10)  into  (8)  gives 

!l  3  (uu)  .  .  b,  3(ulle)  ,  .  h,„  3Ue  ,  ^  3De  , 

-  !  zr  dy  +  1  aT-  dy  +  ,Ve  "SiT  dy  '  1  “S  dy- 

o  o  o  o 

-  and  combining  the  first  two  and  the  second  two  terms  yields 

h  3  h  ... 

/  (uUe  -  uu)dy  +  /  — (Ue  -  u)dy. 
o  o 

Substituting  (7)  and  (11)  into  equation  (6)  gives 

/  (uUe  -  uu)dy  +  /  (Ue  -  u)dy  +  /(Ue  -  u)dy  = 
o  o  o 

The  velocity  gradient  3Ue/3x  is  independent  of  y  as  is  the  partial 

derivative  with  respect  to  x;  thus,  these  can  be  pulled  outside  the 

integration  symbol  so  that 

■1^  /(uUe  -  uu)dy  +  /(Ue  -u)dy  +  /(Ue  -  u)dy  ■  ^  ^ 

o  o  o 

Finally,  rearranging  the  terms,  equation  (12)  becomes 

!r  »«2  f  lb  a-  dy  +  ue  m-  Sr)  ^  +  fr  «•  ^  <13) 


3x 


Ue 


3x 


Ue' 


3t 


Ue7 
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IV  von  Karman-Pohlhausen  Integration 


In  the  preceeding  chapter,  the  theory  for  the  unsteady  flow 
boundary  layer  was  applied  to  develop  the  momentum- integral  boundary 
layer  equation  for  unsteady  flow.  The  steady  flow  version  was. first 
integrated  and  applied  by  K.  Pohlhausen  in  1921.  The  method  as  used 
in  modern  texts,  and  that  is  used  here,  is  a  more  modern  development 
done  in  1967  by  H.  Holstein  and  T.  Bohlen  (Ref.  3:206).  However,  the 
unsteady  flow  application  has  not  been  integrated  because  the  addition 
of  the  transient  term  requires  another  independent  equation  for  closure 
as  will  be  evident  in  the  analysis  that  follows.  The  closure  equation 
can,  however,  be  constructed,  and  the  success  with  this  closure 
equation  is  directly  responsible  for  the  success  of  this  study.  Since 
the  unsteady  flow  method  closely  parallels  the  method  for  steady 
flow,  the  steady  flow  development  of  the  von  Karman-Pohlhausen  approxi¬ 
mate  method  has  been  reproduced  in  Appendix  C. 


The  von  Karman-Pohlhausen  Approximate  Method  for  Unsteady  Flow 
The  momentum-integral  equation  for  unsteady  flow  is 


h Ue2 '  (1-  +  Ue 

o 


m-  fciy  *  If  ue  m-  .  *  <13> 

o  o 


Two  parameters  relating  to  the  boundary  layer  thickness  are  now 


introduced: 


Equation  (14)  represents  the  displacement  thickness,  and  equation  (15) 
represents  the  momentum  thickness  (see  Fig  6).  Substituting  these 
into  equation  (13) : 

3_  Ue2  62  +  Ue  3Ue  .  6L  +  9_(Ue  61)  =  tw 
3x  9x  3t  p 

or  by  expanding  the  first  term 

2  Ue  62  aUe  +  Ue2  362.  +  Ue  3Ue  61  +  3 (Ue  6i)  =  tw  , 

3x  9x  8x  9t  p 

and  finally  rearranging  yields 

(262  +  «i)  Ue  9Ue  +Ue2  362.  +  Ue  Hi  +  61  3Ue  =  tw  (16) 

3x  3x  3t  3t  p 

As  it  stands,  equation  (16)  may  not  be  integrated  in  the  manner 

of  the  von  Karman-Pohlhausen  method.  In  order  to  proceed,  a  closure 

equation  must  be  introduced  which  relates  the  change  in  61  with  time  to 

the  change  in  Ue  with  time.  To  find  the  equation  of  closure,  first 

assume  that  61  is  directly  related  to  the  boundary  layer  thickness: 

«1  =  C16  (17) 

which  is  generally  true  in  the  time  domain.  Then 

361  °  Ci  36 
3t  3t 

As  long  as  the  flow  is  laminar,  it  can  be  shown  that  the  boundary 
layer  thickness,  6,  can  be  related  to  the  velocity  of  the  potential 
flow,  Ue,  by 

6  =  C?/\T~ 

/  Uex 


361  =  ClC2jT  1_  (C^"1/2 
3t  f  x  3t 

=  C1C2  fT  •  1  (Ue|-3/2  3Ue 

Vx  “2  3t 


so  that 


C  c  (-  -i- 

12  Uex  \  2Ue  3t  / 


JllDe  ,  n 

~  ~  2Ue  3t  1 

as  long  as  the  time  dependence  of  and  Cj  are  small  as  compared  to 

that  of  Ue,  which  is  somewhat  restrictive.  Equation  (18)  thus  forms 

the  closure  equation  necessary  to  put  the  unsteady  integral  equation 

in  manageable  form.  Substituting  equation  (18)  into  equation  (16) 


yields 


/nr  .  r  \  3Ue  .  „  ?  862  .  r  3Ue  1  „  3Ue  _  tw  ,,„x 

(2*2  +  V  Ue  +  0e si  +  {1  “at  -  2  S1  It  -?  •  (19) 


Equation  (19)  is  then  multiplied  through  by  62/v: 


(2  +  2<$2.  1«2  +  I /«L\«2l  ^  = 

\  «2/  v  3x  Ue  v  6x  2  V  «2  /  v  3t  y 


and  divided  through  by  Ue: 


( 2  +  |l\  5SS.  +  Ue  &  +  -i- 

l  62/  v  3x  v  3x  2Ue 


1_  (  Sited  25s.  . 

Ue  \«2/v  3t 


T«&2 

yUe 


Finally: 


+  44^ 


_L  lUe 

Ue  3t 


is  added  and  subtracted  to  get 


(2  f  *1  \622(3Ue  +  _1_  3Ue\  Ue6262'  /2  +  I  «l\ «d  X  2HS.  =  3 

V  62  )  v  \  3x  Ue  3 1 /  v  \  2  62  /  v  Ue  3t 


Tws2 

yUe 


To  make  a  boundary  layer  computation,  some  information  must 
first  be  determined  concerning  6^  and  however,  these  parameters 
are  in  terms  of  an  as  yet  undetermined  velocity  profile  within  the 
boundary  layer.  The  velocity  profile  must  satisfy  two  constraints: 
a.  The  order  of  the  velocity  equation  must  be  of  an  order 
compatible  with  the  number  of  boundary  conditions. 
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b.  The  velocity  profile  must  allow  for  an  inflexion  point 


since  the  flow  will  be  with  a  pressure  gradient  (Ref.  8). 


y  =  0  u  =  0  (22a) 

v92u/3y2  =  -  Ue  3Ue/3x  +  3Ue/3t  (22b) 

y  =  5  u  =  Ue  (22c) 

3u/3y  =  0  (22d) 

32u/3y2  =  0  (22e) 


Note  that  except  for  the  second  boundary  condition,  these  are  all 
identical  to  the  steady  flow  boundary  conditions.  The  second  boundary 
condition  is  altered  only  to  account  for  the  unsteady  pressure  term  from 
Euler's  equation  as  developed  in  the  Appendix. 

The  existence  of  five  boundary  conditions  allows  for  a  fourth 
degree  polynomial  expression  for  the  velocity  profile: 

—  -  A  +  Bn  +  Cn2  +  Dn3  +  En4,  (23) 

Ue 

where  n  =  y/6. 

Applying  (22a)  gives 

A  -  0,  (24a) 

and  (22c)  yields 

1  =  B  +  C  +  D  +  E.  (24b) 

Now,  applying  (22d) 

0  =  B  +  2C  +  3D  +  4E  (24c) 

and  applying  (22e) 

0  =  2C  +  6D  +  12E.  (24d) 


Finally,  equation  (22b)  is  applied  to  get 


“  ■ B  +  20,1 +  3D”2  +  ^ 

P  ' 20  +  6D"  +  12E"2|n-0 


_  -62 /3Ue  J_  9Ue\ 
L  v  V  3x  Ue  3t  / 


(24e) 


At  this  point,  a  dimension  less  parameter.  A,  is  introduced  and 


defined  as: 


_  6^/aUe  J_  3Ue\ 

A  ~  v  \  3x  Ue  at/ 


Notice  that  A  here  is  different  from  that  defined  for  steady  flow. 
Substituting  (25)  into  equation  (24d)  and  solving  the  equations  (24a-d) 
yields 

A  =  0 

B  =  2  +  A/6 
C  =  -  A/2 

D  =  -2  +  A/2 
E  -  1  -  A/6 

These  are  exactly  the  results  obtained  for  the  steady  flow 
boundary  layer  velocity  profile.  This  should  not  be  surprising  in 
that  one  should  expect  to  recover  the  steady  flow  solution  from  the 
unsteady  flow  solution  if,  in  fact,  the  unsteady  flow  terms  are  dropped. 
This  is  what  happens  if  the  unsteady  term  is  dropped  from  the  definition 


The  details  of  how  to  solve  for  the  constants  are  in  Appendix  D. 
The  velocity  profile  becomes: 


^  =  ^2n  -  2n3  +  n^)+  A/6  |  n  -  3n2  +  3n3  -  n^) 


and  6^  and  62  can  now  be  solved: 


For  y>6,  u=Ue  by  the  definition  of  boundary  layer,  and  integration 
beyond  that  point  gains  no  additional  information.  Thus,  the  integration 
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to  h>6  gives  a  result  identical  to  an  integration  only  to  h=6.  Thus, 


?! 

h 

i 


n 


R  « 


H  -  -  ujW- 

Substituting  n  into  the  momentum  thickness  and  adjusting  the  integration 
limits  accordingly  yields 

and  substituting  for  u/Ue  gives 

=  /|l-(2n  -  2n3  +  n^)  -  (n  -  3n2  +  3n3  -  dn; 


61  3  A 

6  ~  10  “  120 


(26) 


and  similarly, 

1 
=/ 

0 


^  =f  £(2  -2n3+nS  +  |-(n-3ri2+3n3-ri^)  j*  £l-(2n-2n3+nS  -  £(n-  3n2+3n3-nS  jdn 


62. 37 _ A _ r 

6  315  '  945  '  9072 


(27) 


The  details  of  the  integration  are  presented  in  Appendix  D.  The  shear 
stress  at  the  wall  is  defined  as: 

y3u] 


TW  = 


3y|y=0; 


or 


h£  *  3  (u/Ue)  I  _  7  ,  A 
pUe  9(y/6)jn=0  6 


(28) 


The  separation  criteria  can  now  be  determined,  as  separation  is  defined 
as  the  point  where  the  shear  stress  as  measured  at  the  wall  is  zero  so 
that 


i.e. , 


2  +  A/6  *  0 


A  -  -12 


(29) 
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The  momentum-integral  equation  in  its  last  form  was: 


/ 2+6j\  622  /3Ue+  1  3Ue\  +  Ue6262  -  / 2+1  6l  \  «22  _1  3Ue 
\  62/  v  \  3x  Ue  3 1  /  v  V  2  62  /  v  Ue  3t 


=  tw62  (21) 

yUe 


By  introducing  another  dimensionless  parameter,  K,  defined  by 

K  =  ?./3 Ue  +  1_  3Ue  \ 

\  3x  Ue  3 1  / 


(30a) 


K  ~  6^2  /  3Ue  +  1_  3Ue\ 
v  \  3x  Ue  3t/ 

the  right  hand  side  of  equation  (21)  becomes 


tw62  =  tw6  .  ^2 
yUe  yUe  6 


*  / 

2  +  A\/37  - 

A  - 

A2  \ 

\ 

6/ \315 

945 

9072  / 

and  similarly. 


K  =  6^  (  3Ue  +  1_  3Ue\/ «2.\2  5 
v  \  3x  Ue  3 1 A  6/ 


K-( 


37  -  _A_  -  A2\2  A 

315  945  9072 ) 


Finally, 


(30b) 


|  61  =  61 

.  6  = 

(  ?_  -  JL\/37_ 

-  A  - 

A2  \-l 

(33a) 

62  6 

62 

\  10  120  /\315 

945 

9072/ 

«I  fl(R) 

62 

Substituting  into  equation  (21)  yields 

[2+fi(K)J  K  +  Ue62fi/  -  F2+l,fi(K)j  Z_  3Ue  =  f2(K) 
v  L  2  J  Ue  3t 

Now, 

Ue62fif  =  I  Ue  d(622) 
v  2  v  dx 


(33b) 


Ue62fi/  =  Ue  |dz/dx) 
v  2 


I 
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i 


Substituting  (35)  into  (34)  gives 

r2+f1(K)'l  *K  +  Ue  dz  -  p+lfi(K)]_Z  3Ue  =  f2(K) 

L  2  dx  1.  2  J  Ue  3t 

Ue  dz  -  [4+f^K)]  Z_  aUe  =  2f2K  -  4K  -  2fl(K)  •  K  (36) 

dx  Ue  at 

The  right  hand  side  of  equation  (36)  can  be  defined  as: 

2f  2K  -  4K  -  2Kf i (K)  =  F(K)  (37) 

which  again  is  identical  in  form  to  the  steady  flow  solution.  Indeed 
the  steady  flow  solution  can  be  recovered  if  the  unsteady  terms  are 
dropped.  So  that  substitution  of  equation  (37)  into  equation  (36) 


finally  gives 

dz  = 
dx 

pF(K)  +  [4+f1(K)']Z  3Uel  1 
[_  1  Ue  3tJ  Ue 

(38a) 

Z  = 

K  /  3Ue  +  1  3Ue \-1 

(38b) 

\  3x  Ue  3t  / 

which  are  the  modified  equations  to  be  used  to  solve  the  boundary  layer 
using  the  method  of  von  Karman-Pohlhausen.  Again,  it  is  interesting 
to  note  how  the  steady  flow  solution  is  recovered  if  only  the  unsteady 
terms  are  dropped  from  the  equations. 

As  in  the  steady  flow  solution  method,  the  integration  procedure 
begins  at  the  stagnation  point,  where  Ue=0.  In  this  study,  the  stag¬ 
nation  point  conditions  are  computed  in  the  absence  of  unsteady  flow 
so  that  as  in  the  steady  flow  case  F(K)  must  be  zero  so  that  dz/dx 
at  the  stagnation  point  is  finite.  Solving  (67)  for  F(K)=0,  and  for 
steady  flow  gives 

A0  *  7.052  (39a) 

Kq  -  0.077  (39b) 

Then,  ZQ  -  0.077,  and  dzl  =  0.  Using  the  limit  (Ref.  3:211): 

Uei  dxj^  0 
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Now,  to  complete  the  discussion  on  the  Pohlhausen  integration,  the 


steps  to  integrate  the  boundary  layer  are: 

1.  Compute  the  velocity,  Ue,  velocity  gradient,  9Ue/9x  and 
9Ue/9t  using  the  potential  flow  solution  about  an  airfoil.  See  Appendix 
A  for  details. 

2.  Knowing  zQ  and  dz/dx,  find  z^: 


dz 


Z1  "  Zo  +  dxl*  dx* 


where  dx  is  an  increment  of  arc  length  along  the  airfoil. 

3.  Compute  a  new  value  of  K  from  equation  (38b).  At  this  point, 
the  value  of  A  can  be  backed  out  from  equation  to  determine  flow 
attachment  or  separation  (32).  However,  since  A  for  separation  is  known, 
as  is  K  =  K(A),  then  K  for  separation  can  be  determined  and  thus  allow 
the  elimination  of  one  step  in  the  integration  process.  This  is  done 


for  this  study;  Kgep  =  -  0.1567  (Ref.  3:212). 


4.  F(K)  is  determined  by: 

a)  Knowing  A  and  solving  equation  (32),  (33),  (34) 
and  finally  (37),  or, 

b)  Approximating  F(K)  vs  K  by  a  linear  approximation: 

F(K)  =  a  -  bK 
where  a  =  .47 


b  =  6  (Ref.  3:213) 

Method  (b)  was  chosen  for  this  study. 

5.  Knowing  F(K)  from  the  linear  approximation  and  Ue  and  9Ue/9t 
from  the  potential  flow  solution,  f(K)  by  equation  (33),  and  the  current 
value  of  z,  a  new  dz/dx  is  computed. 


6.  Knowing  z^  and  having  computed  dz/dx|^;  z^+j  is  computed  by 

zi+i  *  z±  +  dz/dxj  ^  .  dx 

7.  Repeat  steps  3-6  with  steady  flow  only  until  the  boundary 
layer  develops  numerically  approximately  a  1-degree  arc  along  the 
airfoil.  Then,  repeat  steps  3-6  using  the  unsteady  methods  until 
the  1/4-chord  is  reached.  Check  K  for  separation.  If  the  flow  is 
attached,  increment  the  angle  of  attack  and  repeat  until  separation 
occurs.  Then  choose  a  different  pitch  rate  and  repeat  again  until 
separation  occurs.  As  separation  is  reached  each  time,  record  velocity, 
pitch  parameter,  and  angle  of  attack  at  flow  separation. 
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Results 


Discussion 

No  theoretical  investigation  would  be  complete  without  a  demon¬ 
stration  of  how  well  that  theory  performed  against  measured  results. 

As  discussed  in  Chapter  II,  there  are  two  experiments  on  record  that 
concern  a  pitching  airfoil  and  pitching  airflow,  and  each  one  shows 
a  positive  correlation  between  this  pitching  and  either  the  delay  in 
stall  angle  of  attack  or  the  enhancement  of  maximum  lift  coefficient. 

To  demonstrate  the  validity  of  the  theory  developed  in  this 
study,  the  modified  von  Karman-Pohlhausen  integral  method  was  programmed 
on  the  AFIT/Harris  500  computer  (see  Appendix  E) .  In  the  unsteady  flow 
Pohlhausen  method,  the  velocity  Ue,  the  time  dependent  velocity  variation 
3Ue/3t,  and  velocity  gradient,  3Ue/3x  are  required.  However,  to  deter¬ 
mine  these,  a  specific  geometry  must  be  defined.  To  do  this,  an  11% 
thick  symmetric  Joukowski  airfoil  (J011)  was  chosen  as  a  model  for 
application.  A  Joukowski  airfoil  was  chosen  because  it  has  its  origins 
in  complex  potential  flow  theory,  and  the  necessary  information  (i.e. 

Ue,  3Ue/3x,  and  3Ue/3t)  is  easily  computed  without  requiring  some  of 
the  other  laborious  techniques  required  to  model  the  flow  over  (real) 
airfoils.  Besides,  one  of  the  goals  is  to  apply  the  theory  to  a 
representative  problem  (and  not  necessarily  an  exact  problem)  to  deter¬ 
mine  how  the  theory  works.  For  the  section  of  the  airfoil  being  used 
in  the  analysis,  it  is  for  all  intents  and  purposes  an  exact  problem 
without  having  to  be  concerned  with  the  trailing  edge  cusp. 

To  approximate  the  required  derivatives,  several  methods  can  be 
used;  namely,  forward  difference,  backward  difference,  and  central 


m 


difference.  It  has  been  noted  that  the  variation  in  accuracy  of 
either  of  these  differences  when  applied  in  the  Pohlhausen  integration, 
is  insignificant  (Ref.  4).  Thus,  a  central  difference  was  chosen  for 
3U/3x: 


3U  =  1^+1  -  Uj-1  , 

3x  d^i+i  “  d6i_i 

and  a  backward  difference  scheme  was  chosen  for  the  time  variation 


—  =  Ui._.-Ulr,i  where  At  =  Ml=l 
3t  At  ’  wnere  Ui-  1 


Results 

The  program  was  run  at  three  different  velocities,  and  at  each 
velocity,  four  different  pitch  rates  were  chosen.  These  were  input 
semi-interactively  into  the  program.  The  output  consisted  of  several 
columns  of  data,  such  as  velocity,  velocity  gradient,  etc.  Of  special 
importance  was  the  column  of  values  of  K.  The  computation  ceased  at 
the  quarter  chord,  at  which  point  K  was  checked.  If  K  was  within  a 
few  precent  of  the  stall  indication  (K  =  -  .1567),  the  angle  and  pitch 
parameter  were  recorded  and  plotted  (see  Fig.  9).  The  16  runs  that 
resulted  in  separation  at  the  quarter  chord  are  plotted  in  Fig.  9. 

The  figure  points  out  a  few  very  important  details: 

1.  The  relationship  between  the  delayed  stall  angle 
and  pitch  rate  is  linear,  as  shown  experimentally. 

2.  There  is  no  change  in  the  effect  with  respect  to 
velocity  of  the  freestream,  when  the  pitch  rate  is 
non-dimensionalized . 

These  results,  while  seemingly  subtle,  are  indeed  profound  in 
that  they  allow  us  to  say  that  a  theoretical  tool  for  prediction  of 
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this  phenomena;  namely,  gust  response  of  an  airfoil,  does  exist. 


Another  result,  shown  in  Fig.  10,  is  the  contribution  the 
transient  term  8Ue/3t  in  the  Euler  equation  of  motion  makes  to  the 
behavior  of  the  flow.  In  Fig.  10,  the  shape  factor  for  both  steady  and 
unsteady  flow  are  plotted  against  chord  location  for  a  typical  run.  No 
in  steady  flow,  only  3Ue/3x  contributes  to  the  shape  factor,  K;  whereas 
in  unsteady  flow,  both  3Ue/3x  and  (^/Ue)‘  3Ue/3t  contribute.  Since  at 
stagnation,  K  *  0.0770,  and  at  separation  K  =  -  0.1567,  it  is  obvious 
that  a  less  negative  value  of  K  results  in  better  behaved  (attached) 
boundary  layer  flow.  As  seen  in  Fig.  10,  this  is  exactly  the  effect 
of  1/Ue  *  3Ue/3t — to  make  the  shape  factor  less  negative.  Thus,  a 
third  observation  can  be  made: 

3.  The  transient  term  (VlJe)  ‘  3Ue/3t,  in  the  unsteady 
Euler  equation  of  motion,  has  a  positive  effect  in 
maintaining  an  attached  boundary  layer  resulting  in  a 
delay  in  the  stall  angle  of  attack. 

Comparison  with  Experiment 

As  discussed  in  Chapter  II,  Kramer  formulated  a  relationship  to 
relate  gust  pitch  rate  to  the  delay  in  stall  angle.  This  was  given  as: 

cftmaxdyn  *  C*max8t  +  0.36  ^  a 

where  a  is  in  radians /second.  Referring  to  Fig.  9,  the  results  of  the 
theoretical  study  show  that 

j-Cn 

astalldyn  «  astallgt  +  16.00  •  — —  . 

If  this  is  converted  to  radians: 

asa„„  ■  +  0.28  a. 


To  compare  this  to  the  Kramer  experiment,  the  relationship  must 
reflect  not  the  difference  in  stall  angle  of  attack,  but  rather  the 
maximum  lift  coefficient.  To  get  the  results  obtained  in  this  study 
into  the  form  that  Kramer  had,  let  it  be  assumed  that  the  slope  Cla 
in  the  linear  region  of  the  CJl-vs-a  curve  be  continued  upward  and  to 
the  right  at  the  same  slope  so  as  to  relate  Aastall  and  ACJlmax..  If 
this  is  assumed,  and  classic  airfoil  theory  is  applied  (i.e.  CHa  = 
2Tr/radian) ,  then  the  results  of  this  study  can  be  expressed  as: 

AC£uiax  =  0.878  Ca  , 

U_ 

00 

where  the  factor  of  \  has  been  included  in  the  slope  as  was  done  in 
Kramer's  experiment  (Ref.  1). 

The  slope  is  considerably  greater  as  determined  by  the  theory. 
However,  it  must  be  noted  that  Kramer  did  his  experiments  on  actual 
wing  surfaces  and  not  true  airfoils.  These  wings  had  aspect  ratios 
on  the  order  of  1.5.  Thus,  if  three  dimensional  flow  corrections 
(Ref.  11:2-6  and  2-7)  are  made  to  the  results  of  this  study,  a  compari- 
tive  result  is: 

AC^max  =  0.343  _Ca  . 

u00 

This  gives  a  remarkably  close  result  when  compared  to  Kramer's 
experiment.  However,  the  actual  aspect  ratio  of  Kramer’s  test  wings 
was  eyeballed,  and  the  actual  values  are  not  specified  in  his  report. 
Thus,  the  corrected  slope  for  the  results  of  this  study  could  fall 
within  a  range  of  0.263  to  0.405,  depending  on  the  value  of  the 
aspect  ratio  chosen. 


VI  Conclusions  and  Recommendations 


Conclusions 

1.  There  exists  a  positive  correlation  between  gust  pitch 
rate  and  the  delay  in  the  stall  angle  of  attack — the  so-called 
dynamic  stall  angle  of  attack. 

2.  Existing  steady  flow  theory  can  be  successfully  modified 
to  mathematically  explain  this  phenomena.  In  conjunction  with  this, 
existing  integral  methods  can  also  be  modified  to  model  this  phenomena. 

3.  For  the  von  Karman-Pohlhausen  integral  method  used  in  this 
study,  the  equation  of  closure  exists  to  successfully  complete  the 
transient  flow  modification,  without  introducing  significant  approxi¬ 
mations  as  to  the  relative  dominance  of  some  terms  vs  others. 

U.  With  only  a  slight  variation  in  slope  of  the  respective  curves 
there  now  exists  both  an  experimental  demonstration  and  a  numerically 
proven  analytical  explanation  of  the  physics  of  gust  flow  response 
of  an  airfoil  that  heretofore  did  not  exist. 

5.  Since  the  potential  flow  acts  in  a  linear  fashion,  the 
pitching  airflow  really  does  not  change  the  dynamics.  Therefore,  it 
must  be  the  pressure  gradient  in  unsteady  flow  as  determined  from 
Euler' 8  equation  that  controls  the  behavior  of  the  stall  dynamics  in 
this  environment. 

Recommendations 

As  in  any  study,  the  end  of  this  effort  does  not  also  mean  the 
end  of  the  investigation.  The  end  of  this  study  should  be  the  founda¬ 
tion  for  many  future  efforts.  Some  of  those  efforts  particularly  high¬ 
lighted  by  this  study  are: 
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1.  In  developing  the  equation  of  closure,  it  was  assued  that 
the  time  dependence  of  the  constants  and  C2  was  small  compared  to 
the  time  dependence  of  Ue.  This  approximation  needs  to  be  studied 
further  and  reported  on. 

2.  In  the  execution  of  the  von  Karman-Pohlhausen  method  for 
unsteady  flow,  it  was  assumed  that  the  behavior  of  ^ x 2  was  not 
affected  by  the  time  dependent  nature  of  the  flow.  This  may  or  may 
not  be  true,  and  must  be  examined  further. 

3.  It  would  be  interesting  to  execute  the  same  problem  using 
a  modern  computational  numerical  method  and  to  compare  the  complexity 
of  development  of  each,  compare  the  efficiency  of  the  computer  program 
for  the  numerical  method  to  the  computer  program  developed  for  this 
study,  and  to  comment  on  regimes  of  applicability. 

4.  The  Deekens  and  Kuebler  experimental  investigation,  the 
pitching  airfoil,  has  yet  to  be  mathematically  modelled.  The  two 
major  problem  areas  that  need  to  be  investigated  before  the  total 
physics  can  be  understood  are: 

a)  derivation  of  a  suitable  governing  equation  that 
adequately  describes  the  boundary  layer  when  the 
airfoil  is  pitched;  and 

b)  the  determination  of  the  appropriate  boundary 
condition  or  conditions  as  applied  to  this 
governing  equation. 

5.  In  trying  to  relate  the  results  of  this  study  to  those  of 

Kramer,  it  was  assumed  that  the  slope  ofCJlvs.  o  between  and 

C^maxst  is  the  same  as  is  for  CJu  in  the  steady  flow  re8ime*  This 

may  or  may  not  be  true,  and  this  relationship  must  be  studied,  either 
analytically  or  experimentally,  and  reported  on. 


6.  In  modeling  the  airflow,  the  trailing  edge  (starting)  vortex 
that  is  shed  when  the  airflow  pattern  changes  about  the  airfoil  (as  is  the 
case  in  this  study)  was  not  included.  These  starting  vortices  are  neces¬ 
sary  to  satisfy  Helmholtz'  Vortex  Theorems  and  are  physically  observed 
(Ref.  12:53-55).  In  this  study,  the  effect  of  these  was  assumed  to  be 
small,  but  this  assumption  may  need  further  study  (Ref.  13:379^390  and 
Ref.  14:376-378). 
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Schematic  of  Kramer  Experiment 
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Figure  7.  Fluid  Volume  Element 
Used  for  the  Continuity  Equation 
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*jCa,  Radians 
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Figure  9.  Effect  of  Pitching  Airflow  on  astall 


Figure  10.  Effect  of  Acceleration  Due  to  Pitching  on 
Boundary  Layer  Shape  Parameter 


b.  Momentum  into  and  out  of  a  Control  Volume  Element 


Figure  12.  Fluid  Element  Along  a  Streamline 
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Throughout  this  study,  many  references  are  made  to  the  inviscid 
flow  velocity,  or  the  velocity  of  the  potential  flow,  and  so  on.  In 
this  appendix,  it  is  shown  how  to  derive  that  velocity  by  using  complex 
potential  flow  theory. 

Complex  Flow  Velocity 

The  complex  velocity  about  a  body  is  the  first  derivative  of  the 
complex  flow  potential  that  describes  the  flow: 

v(z)  =  dF(z)/dz  (Ref.  5)  (Al) 

where  w(z)  =  the  complex  velocity  and 

F(z)  *  the  complex  flow  potential 
Now,  to  obtain  the  expression  for  the  complex  flow  velocity  for 
the  problem  of  this  study,  i.e.  flow  about  an  airfoil,  it  is  necessary 
to  first  find  the  velocity  about  a  cylinder  (that  will  transform 
into  the  desired  airfoil) ,  and  then  transform  that  velocity  via  the 
Joukowski  transformation  into  the  corresponding  velocity  about  the 
airfoil.  The  describing  equation  for  the  complex  potential  flow  about  a 
cylinder  is  a  uniform  stream  at  angle  of  attack  is: 

F(z)  -  Uqq  (z  +  R2/z)  +  i  (T  In  z)/2tt  (A2) 

Then,  the  corresponding  velocity  is: 

w(z)  -  %>  (1  -  R2/z2)  +  i  r/2irz  (A3) 

The  Kutta  condition  that  air  flows  smoothly  off  the  trailing 
edge  (Ref.  6:34)  dictates  the  amount  of  circulation  required: 

r  *  4  Uqq  R  sin  a  (A4) 

for  symmetric  airfoils  (Ref.  5).  Substituting  (A4)  into  (A3): 

w(z)  -  Uqq  (1  -  R2/z2)  +  i  (2  Uqq  R  sin  a) / z  (A5) 
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which  described  the  velocity  w  at  any  point  z(x,iy)  on  a  cylinder  of 
radius  R,  with  its  center  at  the  origin. 

The  advantage  to  using  complex  potential  flow  theory  is  that 
flows  due  to  different  phenomena  can  be  superimposed  to  form  one  flow. 
In  the  problem  of  this  study,  it  is  desired  to  simulate  uniform  flow 
about  a  cylinder  at  some  angle  of  attack  in  two-dimensional  space. 

To  generate  flow  about  a  cylinder , a  doublet  must  be  used.  Flow 
at  angle  of  attack  (that  produces  lift)  requires  circulation,  which  is 
provided  via  the  use  of  a  vortex.  The  describing  functions  of  each  of 
these  in  complex  potential  flow  theory  are: 

Uniform  stream:  Uqq  •  z  ^A6.a) 

Doublet:  Uqq  •  R2/z  (A6.b) 

Vortex:  i  (T  In  z)/2ir  (A6.c) 

(Ref.  7:449) 

Now,  to  describe  the  flow  about  a  corresponding  airfoil,  use: 

|w(g)|  =  |w(z)|  /  |dg/dz|  (A7) 

where 

w(g)  =  flow  about  the  airfoil  (or  body) 
w(z)  =  flow  about  the  cylinder 
dg/dz  =  the  derivative  of  the  function  that  transforms  the 
cylinder  to  the  airfoil. 

The  general  transformation  from  a  cylinder  to  a  flat  plate  (airfoil  of 
zero  thickness)  is: 

g  -  z  +  R2/z  (A8) 

So  that 

dg  ■  1  -  R2/z2 
dz 
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To  generate  an  airfoil  with  thickness,  like  that  which  is  being 
used  for  this  study,  the  center  of  the  circle  must  be  offset  along  the 


real  axis,  into  the  flow  so  as  to  generate  an  airfoil  with  thickness 
with  the  proper  orientation: 

z’  =  z  +  XQ 

where  XQ  =  the  amount  of  offset.  Thus,  the  transformation  is: 

^  =  1  -  R2/z ' 2 
dz 

This  formulation  of  how  to  compute  the  velocity  at  a  point 
C(£,  in)  on  the  airfoil  has  been  programmed  as  a  subroutine  in  the 
computer  program  developed  for  this  study  (see  Appendix  E) . 

Joukowski  Airfoil  Development 

The  theory  developed  in  this  study  is  applied  to  a  Joukowski 
airfoil  to  test  the  validity  of  that  theory.  As  mentioned  earlier, 
a  Joukowski  airfoil  is  used  because  the  potential  flow  about  a 
Joukowski  airfoil  is  easily  described  without  having  to  use  various 
other  complex  transformations  although  they  can  be  used  and  applied. 

The  development  of  a  Joukowski  airfoil  has  its  roots  in  complex 
variable  theory.  It  is  created  by  transforming  a  circular  cylinder 
in  one  plane  through  a  transformation  function  into  the  plane  of  the 
airfoil  (Ref.  12:237  and  Ref.  5).  Although  arguments  against 
Joukowski  airfoils  exist  primarily  because  of  its  cusped  trailing  edge, 
it  is  used  here  because  its  leading  edge  characteristics  and  geometry 
closely  parallel  those  of  several  practical  airfoils. 

To  map  a  Joukowski  airfoil,  either  a  graphical  or  a  mathematical/ 
computer  solution  can  be  used.  In  this  study,  a  computer  solution  is 
chosen.  As  stated,  a  Joukowski  airfoil  is  created  through  a  circular 


cylinder  that  will  map  into  the  desired  airfoil.  In  this  study,  it 
was  desired  to  simulate  a  NACA  1102  airfoil  leading  quarter  chord. 
Through  trial  and  error,  it  was  found  that  an  11%  thick  Joukowski 
airfoil  best  did  this. 

Thus,  to  generate  a  symmetric  airfoil,  a  cylinder  was  offset 
from  the  origin  into  the  flow  and  is  described  by  a  series  of  points: 

z(x,  iy)  =  R  exp  (i  *  0)  +  p,  (A9) 

where  R  is  the  radius  of  the  cylinder,  9  is  an  angle  that  locates  a 
point  (x,  iy)  on  a  cylinder  of  radius  R,  and  p  is  the  distance  that  the 
origin  of  the  cylinder  is  offset  into  the  flow,  p  can  be  either 
totally  real  to  generate  thickness  only  or  totally  imaginary  to  generate 
a  flat  plate  cambered  airfoil  or  complex  to  generate  an  airfoil  with 
both  thickness  and  camber. 

Now,  to  create  the  airfoil,  the  circular  cylinder  was  passed 
through  the  transformation  function 

i>2 

C(C,  in)  =  x(x,  iy)  +  z(x>  -ly-y  ,  (A8) 

where  C(£»  in)  describes  the  points  of  the  airfoil. 

A  plot  of  the  airfoil  is  given  in  Fig.  11. 

Euler's  Equation  of  Motion 

To  develop  Euler's  equation  of  motion  in  unsteady  flow,  use  is 


made  of  a  fluid  element  along  a  streamline  (Fig.  12).  The  sum  of  the 
forces  along  the  streamline  are: 


p  ■  pressure 


dn  =  the  differential  length  normal  to  the  streamline 
ds  =  differential  length  along  the  streamline 
dFs  =  force  on  the  side  of  the  element  in  the  direction 
of  the  streamline. 

To  find  Fs  as  a  function  of  pressure,  look  at  the  sidewall 
geometry  (Fig.  12c).  It  is  evident  from  the  geometry  that: 
dFs  *  pdl  sin  0 

But, 

sin  0  =  3 (dn)  ds  •  1 _ 

3s  dl 

dFs  =  P  9 (dn)ds  ; 

3s 

»  • 

the  sum  of  the  forces  =  -  3 (pdn)  ds  +  p  3 (dn)  ds 

3s  3s 


But, 

3  (pdn)  =  dn  +  p  3  (dn)  ; 

3s  3s  3s 

•  i 

the  sum  of  the  forces  =  -  0£  dnds  (A10) 

3s 

Now,  the  momentum  flux  through  this  fluid  element  is  (Fig.  12b): 

M0  -  Mi  +  3_  (Mcv) 

3t 

which  is  equal  to: 

pdnW  +  3_  (pdnW)  ds  -  pdnW  4-  3_  (pdndsV) 

3s  3t 

For  incompressible  flow, 

3p  =  3p  ~  0) 

3s  3 1 

•  • 

the  mass  flux  is: 


pdnds  (fj  (V2)  +  || ) 

Assuming  dn  and  ds  to  be  independent  of  each  other,  the  mass  flux  is 
equal  to: 

Equating  (A10)  and  (All)  through  the  momentum  principle: 

pdnds  (vf  +  ||)=  -  f*  dnds. 

Dividing  by  pdnds,  equation  (A12)  becomes: 

„3V  .  IV  _  _  1  3£  , 
as  3t  p  3s 


(All) 


(A12) 


(A13) 


the  Euler  equation  of  motion  for  unsteady  flow. 
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Appendix  B 


The  momentum-integral  equation  for  boundary  layer  in  steady  flow 
is  reproduced  here.  Those  not  familiar  with  the  derivation  of  this 
equation  as  it  is  done  here  may  wish  to  review  it  before  reading 
Chapter  III. 

The  momentum- integral  equation  is  developed  by  deriving  the 
continuity  and  momentum  equations  for  a  fluid  element  in  the  boundary 
layer  and  then  integrating  the  momentum  equation  through  the  boundary 
layer. 


M-I  Equation  for  Steady  flow 

The  Continuity  Equation.  In  words,  the  continuity  equation 
states  that  the  net  rate  of  mass  flow  out  of  a  control  volume  is 
equal  to  the  time  rate  of  loss  of  mass  within  that  control  volume 
(Ref.  2:10). 

Referring  to  Fig.  7,  this  gives 


h  n 

incut  =  /pudy  +  3 _ /  /pudyjdx  +  motop 

r>  3x'  O  ' 


rain  =  /pudy 


Notice  that  there  is  no  mass  flow  rate  into  the  control  volume  element 
from  the  body  which  satisfies  the  "no  penetration"  condition. 

The  time  rate  of  loss  of  mass  in  steady  flow  is  zero;  therefore, 
the  continuity  equation  for  steady  flow  becomes: 


S  /(udy)dx  +  motop  *  0 
3x  o 
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Momentum  Principle 

The  momentum  principle  states  that  the  sum  of  the  forces  acting 
on  the  fluid  in  a  control  volume  is  equal  to  the  net  rate  of  transport 
of  momentum  out  of  the  control  volume  (Ref.  2:11). 

Referring  to  Fig.  8,  then  the  x-component  of  the  momentum 
equation  may  be  found  by  finding  the  sum  of  the  forces  in  the  x;- 


direction 


h  h  h 

ZFX  =  -TWdx  +  / pdy  -  / pdy  -  3_(/pdy)dx, 
o  o  3xo 


and  equating  it  to  the  momentum  terms : 


h  h  h 

/puudy  +  3_(/puudy)dx  -  /puudy  +  motop  .  Ue, 
o  3xo  o 


such  that 


h  ^  h 

3  (/puudy)dy  +  motop  •  Ue  -  -twdx  -  3_(/pdv)dx.  (B2) 

3x  o  3x  o 

m0top  can  be  replaced  with  its  solution  obtained  in  the  develop¬ 
ment  of  the  continuity  equation 


“otop  “  -3_(/pudy)dx,  (Bl) 

v  3xo 

so  that  the  momentum  equation  becomes 

h  h  h 

3  (Zpuudy)dx  -  Ue  .  3  (/pudy)dx  ■  -twdx  -  3_(/pdy)dx. 

3xo  3xo  3xo 

All  the  integrations  are  with  respect  to  y.  Since  x  and  y  are 

Independent  functions,  the  partial  derivatives  can  be  pulled  inside  the 

integrals: 


h  h  h 

/  3  (puu)dydx  -  Ue  .  /  3_(pu)dydx  ■  ~rwdx  -  /  3  (p)dydx 
o  3x  o  3x  o  3x 


From  the  boundary  layer  assumptions,  it  is  known  that  3p/3y  is 


zero.  Since  this  is  true,  and  since  the  problem  is  two-dimensional,  the 
pressure  gradient  in  the  x-direction  is  the  only  pressure  gradient,  and 
the  partial  derivative  of  pressure  can  be  replaced  by  a  total  (not 
substantial)  derivative: 

h  .  ^  h  . 

/  (puu)dydx  -  Ue  .  /  —  (pu)dydx  =  -  rwdx  -  /  ^  dydx  (B3) 

o  o  o 

If  the  fluid  is  further  restricted  to  incompressible  flow, 
then  p  can  be  divided  through.  Also,  divide  through  by  the  common  dx 
term: 


/  -~(uu)dy  -  Ue  .  /  ~<u)dy  =  -  ?*  -  ±  f  £ 


Now,  Euler's  equation  for  steady  flow  can  be  used  to  substitute 
into  equation  (B3)  for  the  pressure  gradient: 

_  I  =  Ue  (B4) 

p  dx  dx 


h  A  h  _  h  jf. 

3  /  \  ,  3  •  Tu  ...  dUe 


/  — (uu)dy  -  Ue  .  /  —  udy  =  -  +  /Ue 

dX  dX  p 

O  O  O 


dx 


dy 


Rearranging  terms: 
h 


,  Juu  .  .  „  ^  3u  .  .  dUe 
-  /  —  dy  +  Ue  /  ^  dy  +  /Ue 

O  O  O 


(B5) 


Notice  that  the  second  term  looks  like  part  of 

h 


3 (uUe)  dy. 


3x 


which  is  a  form  that  is  identical  to  the  first  term.  Expanding  this 
term  out 


3 (uUe) 


3Ue 


3u 


dy  “  / u  -r-  dy  +  /Ue  —  Cy.  (B6) 


3v 


The  second  term  on  the  right  hand  side  is  the  second  term  in 
equation  (B5).  Thus,  if  the  first  term  on  the  right  hand  side  of 
equation  (B6)  is  both  added  and  subtracted  from  equation  (B5)  (hence 
adding  zero  to  the  equation) ,  it  looks  like 


~f  +  Ue  .  /  dy  +  /u  dy  +  /Ue 


dUe 

dx 


h 

dy  -  fu 
o 


dUe 

dx 


dy 


lit  (B7) 

P 


Combining  the  second  and  third  terms  of  equation  (B7)  to  get  the  term 
of  the  left  hand  side  of  equation  (B6)  gives 


- 1  fe(uu)dy  +  1  fe(uUe)dy  +  'Ue  ^  dy  - /u  dy "  T  (B8) 

o  o  o  o 

I 

The  velocity  gradient  at  the  edge  of  the  boundary  layer.  dUe/dx 
is  a  function  of  x,  along  the  direction  of  flow,  and  is  independent  of 
the  height  away  from  the  body.  Thus,  this  term  can  be  brought  outside 
the  integral.  In  addition,  the  first  two  terms  and  the  second  two  terms 
are  similar  and  can  now  be  combined: 


/  |^(uUe  ~  uu)dy  +  /  (Ue  -  u)dy  * 

o  o 

Bring  the  partial  derivative  outside  the  integral  and  rearrange  the 
terms  to  get 


h Dc2 '  k  l1  -  feb +  is  '(L  -  feW  ■ 1 


1st 

p 


(B9) 


which  is  the  momentum-integral  equation  for  boundary  layers  in  steady 


flow. 


Appendix  C 


The  von  Karman-Pohlhausen  integral  method  for  steady  flow  is 
reproduced  here  again  for  those  readers  not  familiar  with  its  development. 


Method  of  von  Karman  and  Pohlhausen  -  Steady  State 

The  momentum- integral  equation  for  steady  flow  is: 


(1 

3  Ue2  // u_ \  /  _  u_\dy  + 

3x  o\Ue/\  "  U e/ 


Ue  dUe  /  /.  u  \dy  =  Tw  (B9)  and 

dx  o  \  Ue/  p  (Cl) 


Two  parameters  relating  to  the  boundary  layer  thickness  are  now 


introduced: 


h 

«1  »  / 


!  (l  _  «L.\ 

o  \  Ue/ 


«2  -/W_  n  iLNdy  (c3^ 

o  Ue  '  Ue' 

Equation  (C2)  represents  the  displacement  thickness.;  equation 
(C3)  represents  the  momentum  thickness  (see  Fig.  4).  Substituting  into 


the  momentum- integral  equation: 


3  (Ue2fi2)  +  Ue6i  dUe  =  tw 
3x  dx  p 


Expanding  out  the  first  term: 


2Ue52  dUe  +  Ue2  d62.  +  Uedj  dUe  -  TW 
dx  dx  dx  p 


(262  +  61)  Ue  dUe  +  Ue2  d62  *  In 
dx  dx  p 

For  the  boundary  layer  to  be  computed,  then  something  must 
be  known  about  61  and  62*  However,  these  are  in  terms  of  an  as  yet 
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unspecified  velocity  profile,  u/Ue.  The  velocity  profile  must  satisfy 
two  constraints: 


* 

•  ,* 

* 

i 


e 


1.  The  order  of  the  velocity  equation  must  be  in  accord  with 
the  number  of  boundary  conditions. 

2.  The  velocity  profile  must  allow  for  an  inflexion  point, 

since  the  flow  will  be  with  a  pressure  gradient  (Ref.  3:207).  . 

The  boundary  conditions  available  are: 


y  =  0  u  =  0  (C6.a) 

v32u/3y2  =  -  Ue  dUe/dx  (C6.b) 

y  *  6  u  *  Ue  (C6.c) 

3u/3y  =  0  (C6.d) 

32u/3y2  =  0  (C6.e) 


The  existence  of  five  boundary  conditions  allows  for  the 
solution  of  five  constants.  In  other  words,  a  4th  order  polynomial  is 
required  for  the  velocity  profile.  The  form  of  the  velocity  profile  is: 

-  A  +  Bn  +  Cn2  +  Dn3  +  En^ 


Applying  (C6 . a) : 


Applying  (C6.c): 


Applying  (C6.d): 


Applying  (C6.e): 


A  -  0 


1  =  B  +  C  +  D  +  E 


0  -  B  +  2C  +  3D  +  4E 


0  =  2C  +  6D  +  12E 


(C7.a) 


(C7.b) 


(C7.c) 


(C7.d) 


Applying  (C6.b): 


36(u/Ue) 


B  +  2Cn  +  3Dn 


2 


+  4En3 


.  (u/Ue)  =  2C  +  6Dn  +  12En2l 
3'(y/6)2  n  =  0 

2C  =  dUe  (C7.e 

v  dx 

At  this  point,  a  dimensionless  parameter,-^-,  is  introduced: 

A  5  dUe  (C8) 

v  dx 


Substituting  (C8)  into  (C7.e)  and  subsequently  solving  equations  (C7.a) 
through  (C7.d): 

A  =  0 


B  =  2  +  A/6 
C  =  -  A/2 
D  -  -2  +  A/2 
E  -  1  -  A/6 


For  details  on  this  solution,  see  Appendix  D. 
The  velocity  profile  is  now: 


u  *  (2n  -  2n3  +  ri^)  +  A/6  (n  -  3n2  +  3n3  -  n*) 
Ue 


and  Aj  and  $2  can  now  b®  solved: 


%  ’  i  HI  -  Sr)dy 


6  6  *0'*  Ue' 


For  y>  6,  u  ■  Ue  by  definition  of  the  boundary  layer,  and  integration 
beyond  y  *  6  makes  no  sense.  Therefore: 


■  i  i  «  -  6s>dy 


Previously,  n  was  defined  as  y/6. 


Introduce  this  into  the 


displacement  thickness,  and  redefine  the  limits  of  integration 


Substituting  for  u/Ue,  this  gives: 


«r 


— g  =  1  -  (2n  -  2n3  +  n^)  +  A/6  (n  -  3n2  +  3q3  -  n^)Jdn; 


Integrating  this  out  and  evaluating  at  the  limits  yields 

il  -  3  A 

6  “  10  120 

Similarly,  for  the  momentum  thickness, 

1 


(C9) 


02 


=  f  ^(2n  -  2rj 3  +  nS  +  (n  -  3n2  +  3n3  -  n^)^|. 

[l  -  (2ti  -  2n3  +  n*)  £  (n  -  3n2  +  3n3  -  dp. 


Integrating  this  out  and  evaluating  this  at  the  limits 

Al  JL  JL  A2 

6  *  315  945  “  9072  * 

For  details  of  the  integration,  see  Appendix  D. 

The  shear  stress  at  the  wall,  t„,  can  be  written  as: 

w 

Ue  3 (u/Ue) 


(CIO) 


w  *  V  6  3(y/«) 

*  2  +  A/6. 


tb  m  3 (u/Ue) 
yUe  *  3(y/6) 


Cll) 


This  also  gives  the  separation  criteria,  since  separation  is  defined 
as  when  the  shear  stress  at  the  wall  becomes  zero.  Applying  this 
definition  gives 


^  -  2  +  £  -  0. 

yUe  6 


A  *  -12, 

the  criteria  for  flow  separation. 

Rearrange  the  momentum-integral  equation  into  a  more  workable 

form: 


(C12) 


(2S2  +  42|  +  Ue2  ^ 


(C5) 
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Multiply  by  62  and  divide  by  v: 


6?2  (2  +  61/62)  Ue  dUe  +  Ue2  6262.'=  tw52 
v  dx  v  \> 


divide  by  Ue: 


(2  +  02—  iMS.  +  Ue  6262*=  tw62 

S2  v  dx  v  Uev 


Another  dimensionless  parameter  is  now  introduced: 


(C13) 


K  =  z  dUe/dx 
K  =  622  dUe/dx 


(CIA. a) 
(CIA .b) 


In  equation  (C13),  the  right  hand  side  is: 

Twd2  =  TyS  .  62 
pUe  pUe  6 

which  can  be  expressed  strictly  as  a  function  of  as 


tw52  - 


*«>(& 


37  -  A  -  A2 

315  9A5  9072 


which,  for  the  Pohlhausen  technique,  is  defined  as  f2(K). 
Similarly: 

K  =  A  (§2.\2 


(C15.a) 


Finally: 


=  /  37  -  A  -  A2  \2A 


AIL 

\315 


9A5  9072) 


(C15.b) 


«1  *  il  •  i_  -  (±  -  A  W  37  -  _A_  -  a2  yi  ,rl,v 

62  6  62  U0  120/  \315  120  9072/  1  7 


which,  again,  for  the  Pohlhausen  integration  is  defined  as  fi(K).  Thus, 


the  momentum- integral  equation  becomes: 


[2  +  fi(K)]  K  +  Ue6262* =  f2(K) 


(C17) 


Ue6262' -  1  Ue  d(622) 
v  2  v  dx 
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(C18) 


Ue62$2  =  ^  Ue  d (62^/v) 

v  2  dx 

Ued2^2  =  1^  Ue  dz 

v  2  dx 

Substituting  (C18)  into  (C17)  gives 

[2  +  fi(K)]  K  +  1  Ue  dz  -  f 2 (K) , 

2  dx 

or  Ue  dz  =  2f2(K)  -  2Kf1(K)  -  4K  (C19) 

dx 

Now,  the  right  hand  side  of  equation  (C19)  can  be  defined  as 

2f2(K)  -  2Kf 1 (K)  -  4K  =  F(K) ,  (C20) 

so  substitution  of  (C20)  into  (C19)  yields 

dz  =  F(K)  ‘  (C21.a) 

dx  Ue 

z  =  K/Ue'  (C14.a  and  C21.b) 

which  are  the  only  two  equations  required  to  solve  the  boundary  layer 
by  the  method  of  von  Karman  and  Pohlhausen 

The  integration  procedure  begins  at  the  stagnation  point  of 
the  airfoil.  Recall  that  at  the  stagnation  condition,  Ue  *  0,  to 
maintain  a  finite  dz/dx,  F(K)  =  0.  Solving  equation  (C20)  with  appro¬ 
priate  substitutions  for  fi(K),  f2(K.)»  and  K,  it  is  found  that: 


Ao  »  7.052 

(C22.a) 

Ko  *  0.0770 

(C22.b) 

For  details  on  this, 

see  Appendix  D. 

Thus, 

Zq  «  0.0770/U'o 

and 

dz/dx  ^  =  0 

Using  the  limit  (Ref.  3:211): 


dz/dx^  =  -0.0652  U’’0/(U’0)2 

Now,  to  complete  the  discussion  on  the  Pohlhausen  integration, 
the  steps  to  integrate  the  boundary  layer  are: 

1.  Determine  the  velocity  and  velocity  gradient  dUe/dz  from 
the  potential  flow  solution  for  the  shape  of  interest  (see  Appendix  A) 
in  addition  to  the  arc  length  segments  along  the  body. 

2.  Knowing  z0  and  dz/dx, 

zi  =  z0  +  dzl .  dx, 
dxo 

where  dx  is  an  increment  of  arc  length  along  the  airfoil. 

3.  Compute  a  new  value  of  K  from  equation  (C21.b).  At  this 
point,  the  value  of  A  can  be  backed  out  from  equation  (Cl5.b)  to 
determine  flow  attachment  or  separation;  however,  an  easier  way  to 
determine  flow  separation  (and  the  method  used  for  this  study)  is  to 
find  the  value  of  K  that  corresponds  to  A  =  -12,  and  use  this  as  the 
flow  separation  criteria.  The  value  is  K  =  -1567  (Ref.  3:212). 

4.  F(K)  is  determined  either  by: 

a.  Knowing  A  and  solving  equation  (C19) ,  or 

b.  Approximate  F(K)  vs  K  by  a  linear  approximation: 

F(K)  -  a  -  bK  (C23) 

with  a  =  . 47 

b  =  6  (Ref.  3:213) 

5.  Knowing  F(K)  from  the  linear  approximation  and  Ue  from  the 
potential  flow  solution,  determine  the  new  value  of  dz/dx. 

6.  Compute  a  new  z  by: 

zj.+l  =  zi  +  d_z  .  dx£ 
dx  i 

7.  For  this  study,  repeat  steps  3  through  6  until: 
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a)  The  1/4-chord  point  is  reached.  Check  K  against 

the  separation  criteria.  If  the  flow  is  not  separated 
(K  =  -0.1567),  select  a  higher  angle  of  attack,  and 
repeat  until  flow  separation  at  the  1/4-chord  is 
indicated. 

b)  Then  stop  the  computations  and  record  the  steady 


state  stall  angle. 


Appendix  D 


Solution  of  the  Constants  in  the  Velocity  Equation 

The  three  equations  that  remain  to  be  solved  are: 

B  +  2C  +  3D  +  4E  =0 
2C  +  6D  +  12E  =  0 
B+C+D+E=l. 

Substituting  C  =  -A/2,  and  setting  up  an  augmented  matrix 

“  1  3  4  |  0  A 

0  6  12  |  0  A 

_  1  1  1  |  1  A/2_ 

Now  performing  elementary  row  operations: 

Row  2  by  6: 


1 

3 

A  1 

0 

A 

0 

1 

2  1 

0 

A/6 

1 

1 

1  ! 

1 

A/2 

Row  3  -  row  1  +  2  •  row  2 : 

fl  3  4  I  0  A 


L0  1  2  |  0  A/6 

0  0  1  |  1  -A/6_ 

Row  1  -  4  ■  row  3;  row  2  -  2  •  row  3: 

~ 1  3  0  |  -4  +  5/3  •  A 

0  1  0  |  -2  +  A/2 

0  0  1  |  1  -  A/6 


Row  1  -  3  •  row  2: 


1  0  0  |  2  +  A/6 

0  1  0  |  -2  +  A/2 

0  0  1  |  1  -  A/6 
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Thus,  the  constants  are: 


A  =  0 

B  =  2  +  A/6 
C  =  -  A/2 

D  =  -2  +  A/2 
E  =  1  -  A/6 

Solution  of  6i/6: 

^  '  £(l  -  Bi)d" 

*  ^C1  ■  (2 + t)n  +  (■  f)  "2  +  (_2 + ?  y  *  (2  -  hy  ]“*" 


\n2  _ 

(-  AP3  - 

1-2  +  At!4 

It  A\n5 

n 

l  2  Ji 

r2  +  2  h 

l1  6)5 

■(1-1  +  ?-5)+A(-n  +  i-?  +  ^ 

6j_ _ 3 _ A_ 

6  10  "  120 


Solution  of  62/6: 


Az 

6 


u 

Ue 


i'-k) 


dn. 


It  is  more  convenient  to  solve  62/^  in  terms  of  the  constants  A,  B,  C, 
D,  and  G,  and  then  substitute  for  their  values  later. 


§1 

6 


/Fbh  +  Cn2  +  Dn3  +  En^Hl  -  Bn  -  Cn2  -  Dn3  -  En^)~|dn 
0U  J 


/fBn  +  (-B2  +  c)n2  +  (-2BC  +  D)n3  +  (-2BC  -  c2  +  E)n4ldn 
0L  J 
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to 


=  [(-2BE  -  2CD)n5  +  (-2CE  -  D2)n6  +  (-2DE)n7  +  (-E2)n8^]dn 

Bn2  ,  (-B2  +  C)  ■*  .  (-2BC  +  D)  l  .  (-2BD  -  C2  +  E)  Q 
=  -y-  +  - ^ - 1  +  - 4 - n  +  - 5 - h3 

+(-2BE  -  2CD)r6  +(-2CE-  D^,  +  ^8  +  Z§V  1 
6  7  8  9  Q 

Substituting: 

B  *  2  +  A/6 
C  =  -  A/2 

D  **  -2  +  A/2 
E  =  1  -  A/6 

=  .1175  -  . 001A  -  .OOOllA2 

0  ’ 

=  37 _ A_  A2 

315  “  945  "  9072 


Solution  of  Stagnation  Condition  for  A. 

At  the  stagnation  point,  F(K)  *  0  so  that  dz/dx  at  that  point 
is  finite,  since 


(Ref.  3:210) 

Setting  this  to  zero: 

(&  -  m  -  5572)[2  -  Ml A  +  mY2  +  9A2  a3]  ’  °-  (D1) 

This  can  be  expanded  out,  and  solved  by  mathematical  tools, 
or  we  can  try  different  values  of  A  until  equation  (1)  is  satisfied. 

This  gives: 

Ac  -  7.052; 

and 

Kq  -  0.0770. 
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ON  The  C0',b’jT4TICNS  0,'i*NSTSE4M  of  The  stag¬ 
nation  POINT.  ONCE  Twe  BOUNDARY  L  4 x F ,)  HAS 
STE40IE0  OU  T ,  CCMP(J  T  4  T  I  ON  »AS3E3  0  “i  TO  tug 
PITCHING  4I9FL0W. 

NOT  IM’^EO  I  *  TElT  SUBJECTING  Th£  rtfj'  on.llfr  It  fro 
TO  The  PITCHING  AI«FLO*  “AS  JO  JOASf  r  M'li.P 
AFFECT  QN  ‘  "  CONFUTATIONS,  A  S  f‘*F  aoii»|OARY 
LA  yep  IS  ally  YEPY  r  ,t  r  n  fAS  C’^tsf^  t,. 

OOyiNSTPFI  .no  The  GRO-yTm  r.RADIFMT  F«F‘.  j 
TmE  UNSTEOY  FLO  *  IS  ABOUT  T-f  SA‘E  »S  TmF 
GPOajTh  GRADIENT  FO  -  3  T£  40 x  Ego-v,  I*'  r-p  oer,  fo*. 
VERY  NEAP  The  STAG  . ..  T  I  ON  3OI*)T. 

N  a  N  ♦  I 

CONFUTE  the  PEPTI'iEtT  BOUNDARY  LA.fp  j.\o  a  e  ms 

ZZ  *  OZOSE^SSt  ZZ 
RK  *  ZZoOODS 
FK  s  .<17  -  *>,•  RK 
OZOS  a  FK/'Jl 

OELT  a  OSS/Ul 
TINE  a  TI-6  ♦  CELT 
OALPHi  a  OELT  •  400 T 
ANGLE  a  ANGLE  ♦  OALPHA 
AL°H|  a  4LPHJ  ♦  OALPha 

CALL  U  (  ANGLE  #  P  A  0 1  US ,  C  ON ,  e  I » •<  I  *<F ,  alRhA..,2) 

OUOT  a  u2  •  U1 

OuOT  a  OUCT/OELT 

ANGLE  a  ANGLE  •  0.9111 

ANGLE!  *  ANGLE  -  9.J9H 

CALL  IJ(  ANGLE  I  »  P  AO  IUS#  CON,EI.:JtNF,  ANU,  *L°N4,  .2) 
CALL  OSCANGLFUPAOIuS.COn.Anl,  <2, r?) 

ANGLEO  a  ANGLE  ♦  0.0101 

'  CALL  »J(ANr,LEO»PAO|US,CO\,EI,  It'  F.AH.,,  il.5-4,'  ->) 

■call  csi  angle  o*p»o  us  .con  »  a«*l,  xo,  <d 

CALL  U{  ANGLE.PAOI'jS.C-N.FI.uI'.F,  a  -n,  alp-a,. m 
CALL  OS  (ANGLE.  PAN  US,  C  On,  a -J,  *t ,  r  |  J 
OS!  a  SOPTim  -  xl)t*2  ♦  m  .  Y1)**2) 

032  «  SQPT<{«2  -  x  1 )  «  *2  ♦  (Y*  -  v  I ) • *2 ) 
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uuuuuuuuuuuuuuuuuuuuu  uuu  uuuu  uu  u 


css  s  ( '‘S i  ♦  -<5p i  /r  *■>- ' 

z  ("2  -  •"  )  /~3  > 

i.Ti'CS  *  I  •  '•  'S 

tr,r  s  (*2  ♦  *LE  ) /C*--  R  ' 

if(n.lt.sooo)  -,0  ro  in 

N  s  o 

*«MTF(fc,n*0C.  it,  JM.'S,  '"'DT,  F<,  /•->; 

10  CONTINUE 

c 

AOnT  a  4C0T\ 

N  a  0 

00  20  J  *  <l,*2 

THIS  coop  COMPUTES  ThF  /  [no  )(•'  TjO  ^  .  :  -•  y 

LAYEB  »s  IT  !S  SUBJECTED  TO  4  0( r C - 1 '■  5  iNfm  . 

THE  TRANSIENT  FLO*  DYlA-ICS  A-F  O'3  T  oo- 
I*iCqEM6\ T 1 00  T^E  FLOn  A'CLE  *  T  ESC-*  ="I'T  -s 
C0MPU  T  4  T  1 0% ,  t,  BY  4',  AMOUNT 

OALPHA  a  (PITCH  RATE)  •  (T['-P  J'.rBTvr  T) 

AT  EACH  POINT, THE  BOUNDARY  LArE*  pjvvfteaS 
APE  COMPUTED,  AS  APE  r^E  ARC  LF'GT-,  local 
velocity,  TI“E  a vo  SPATIAL  velocity  de¬ 
rivatives,  ANQ  T I ve  tvCRE^AT,  T**E  > 

are  carried  out  from  just  past  t-*e  s t  a g*  a  r  t 
POINT  UP  TO  The  QUARTER  C“*0R0  »0[NT  on  r><t 
SUBPACE  OF  THE  AIRFOIL,  AT  WHICH  301*  f  ThC 
LOOP  IS  EXITED.  r'HpjtATIOV  TO  A  CIFt-ctfT 
POINT  CAN  SE  AC  I  3Y  'PPFLY  C-A-Ol'G 

THE  TEST  CONOI1 4 

N  s  N  ♦  t 

COMPUTE  the  PERTINENT  BOUNDARY  LATpR  baravctprb 

ZZ  a  OZOS«DSS+  ZZ 
RK  3  ZZ»(0UDS  ♦OflDT/Ul  ) 

FK  a  0 ,  a  7  -  6,*SK 
CALL  POHL(RK,RLAMOA J 

DEL2  a  37. /SIS.  -  SL»<OA/RtS.  -  {  *<L A* DA  * • i t/41 72  . 
F2K  s  (.J  -  RLA«nAm.)/CEL2 
OZDS  a  (FK  ♦  (a  ♦  F  2< )  •  ZZ  •  Oii)T/Ll)/i  1 
OZOS  a  FK/Ui 

COMPUTE  The  TIME  INCRE“ENT  fob  a  paoticlF 
TO  TRAVEL  from  POINT  (<)  Tn  POf*iT(l*t) 

CELT  a  OSS/Ul 
TIME  a  time  ♦  OELT 
OALPHA  a  OELT*  AOOT 
ANGLE  a  ANGLE  ♦  OALPhA 
,  ALPmI  a  AL°Hl  ♦  DALPHA 
CALL  U(  ANGLE  ,BAOIliS,CON,EI,UI,f,A'*",ALPhA,.  ri 

COMPUTE  THE  UNSTEADY  VELOCITY  GRADIENT 

OUOT  a  (U2  -  UD/DELT 
ANGLE  «  ANGLE  -  0.0» 

ANGLE  1  a  ANGLE  -  0.0 1 

CALL  U(A»iGLEl,RA*r.S,CP*',Fi,u!\F,  a>«.',al 2) 

CALL  OS(A*iGLEt  ,»AntuS,CDN,A«*.j,*2,rR) 

ANGLED  a  angle  ♦  0.1| 

CALL  Of  ANGLER, RA^ri. k.Cdn.E  1,'V  F,  »•  ,  mr.-a,  -1 

^11  I  7  A  t  P|  C  S  JA'ti  1  P-'»  Y  *  Y  »  ' 


w._ 


nono  n  r>  r>  o  noo 


CM.L  A  -L-  1,  -  4  '  r ■  }.r  ■  ,  •  ,  .  , 

c»tu  vr 4'i',L*»;a  't  s.c  t  * .  ■  .  m. 

Caul.  r.(.r,RAr  ,v>, 

Cp.MPnrF  »or  l£‘-GT^  a*.  ,  ,Fl'CTT»  (r  •  f 

OS l  a  S09T((<1  -  ♦  (ii  -  i 

0S2  *  S09T((<2  -  <1)»*2  ♦  ry?  - 

OSS  a  (OS1  ♦  OS2W0-O5O 

O'jOS  a  (U2  -  >„0)/0SS 

’.’DUOS  *  VI  *  0"C3 

*0C  a  X2  ♦  XLE 

XOC  *  X0C/CHO9O 

STOP  TH£  COMPUTATION  4T  Tip  <3ii4STE’-C'-f'l:o 

IFfXOC.GE.0.2S0)  GO  TO  2S 
IP<*.I.T,2S0)  GO  TO  20 
MaO 


20 

*9 1  Tf  ( 6,  1 )  XOC  * 
CONTINUE 

Ul, 

O'JOS, 

C»)nT. 

,  o<,  ll. 

f^rs 

25 

-pitfos.  n*oc. 

01, 

p'ns. 

9'jor, 

9P.  77, 

'/"S 

**«ITE(6,U5)4LPH1 
*»I TF(6,S5)PtTCH 
*9ITE(6.60)9« 

"prTF(6,7<no€i. 

«»ITE(6,0O) 

Vi9ITE(6,8nTI«E 

1  FP9-4T{«X,P6.  J,2C«*.FtO.SJ,ax,F7.S,a*,F;.  ),?(  j<(-  M)  1 

SO  E09'*»T(tHl,"90t'NfH9r«L*»E9  ?  vj  4 -p  T  r 3  c~o  •<  #  ?*,  #  ^  f  "c  r  /  -\e  <•  « 1  \ 

«0  FOP^tTC"  INJTjai.  ANGLE  CF  ATT'P,.  •’^‘>,3,"  e;-ecs*/) 

45  FC9V»T(/*  final  angle  pe  1  r  r  1  c  1 :  ••-g'>sp«" /1 

SO  E03''AT(»  PITCH  9ATE;  •'PGiFfcS/S£CH/) 

SS  FpOH4  f ( *  PITCH  04H4«ETfcH5  %F7.S/> 

60  F09**A  T  (  •  X  AT  THF  GU  H  TEP-C-  jnn  !  "»F4,<J/l 

i0  F 09MAT(*  fl«L  THIC<N€3S»  FT"/) 

80  EO»**AT("  T  T  VE  TO  9E4CH  THE  GiJ4»  TF^.«  3 

81  FC9-A  T  ( •  CH09Q  F90*«  ThF  S  T4G*’4  T \f>-,  a-<I‘.Ts  »,F7.5,“  ?cCh/> 

STOP 

EMO 


* 


<*3> 


^  K°  ”  D" 
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nnonn  nrinnnnnnnnnnnorinooooo 


SUa»0t'Tl^E  '■!(  ANGLE. -MOT  S,C*- '  e.  >■  '• .  ' 

CO«Pt EX  C«PLX,?,FI, jrrt<,  • 

THIS  SUBROUTINE  !S  USED  T)  zr"  pi.  r£  r-.r  f'fL 
value  OF  T^E  VELOCITY  0*.  i  jr  «r:.c  ...  . 

USING  COMPLEX  °OTC'irtfti.  F'_o«  t'-c«!}v,  »•  T  •  ’  > 

THF09Y,  T >-E  VELOCITY  IS  ff’ST  rf'iO"tri'  js*-  T 

A  CIRCLE  IV  THE  C^PUIATI  =  LA*E;  T-f 
EMVING  EQUATION  IS  The  Of  -  I  .  4  T  T  VF  ->F  T-F 
COMPLEX  FLO«  POTENTIAL  EGrATIO'.  G*  CF  W* 

VELOCITY  has  BEEN  0£TEa«I'E'',  T-c  '.fl'CITy 
ABOUT  The  CCROESP0N0PG  AIBfiil  t3  C'v’°<.  r'-' 

BY  | 

U(AISF0IL1  *  U(CIRCLF)/T ‘?/r;erA) 

xheBE  DZ/OZETa  IS  The  CE«IVATTvE  ’F  T-T 
equation  useo  to  t9amsf*c«  fp.-.v  t~c  r i *•* c i ^ 

TO  The  4I9F0IL. 

the  COORDINATES  »N0  VELOCITY  Im  T-F  Cl^ci ..  f 
PLANE  A9E  COv"UTEO> 

X  a  RADIUS  *  COS(ANGLE«CQ  .) 

Y  a  RADIUS  •  SIN(A\'C-LE  *C0'O 
Z  a  C“PLX(X,Y) 

W  a  1.1 1 *IF  »((!,, 0.)  -  (RADIUS**?) /T«*?  ♦ 

I  (2,*  El  •  RADIUS  *  S  I  v  (  ji.PwAtCovii  )/Z) 

X  a  X  t  AMij 

Z  IS  nOa  rHA»GEO  TO  REpoFSr’  T  T“E  al'-fs 
OF  THE  COOP'Oi*iATES  USED  IV  The  Tna«  s- 
FOShatION  EQUATION* 

Z  a  C«pL*(X,Y) 

OZETA  a  (Z»*2  -  R40t'JS»*2) /Z**2 
UU  a  C4B3(hJ/C4BS(0ZEtA) 

RETURN 

FNO 
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SUBROUTINE  DSC  *Nr,L£,  940tiJS.C^'i,«  -  ,  t .  O 
COMPLEX  C-Pl*»Z 

THJ3  SUBROUTINE  TS  USED  TO  Cr-P’.TE  T -E 
»»C  LENGTH  (>ETv»EEM  r«0  SPECIFIC  l s*  T  S 
ON  T«g  JOL'KO^SKI  4IPFOIL.  3E  t?  -i"1?  ~F 
COMPLEX  V  49 1  48LE  THEORY  TO  ryi  Tuts. 

FIRST,  POINTS  OP  INTEREST  4Cf  CC  'P.TPO  r\ 

4  C IPCLE  •  4H0  49E  T  HF  N  TRaNSFOR'-''"'  InT'' 
C0090IN4TES  OP  THE  4  I PF  0 IL  5  V  Tup  P0i'4Tt<-": 

ZET*  S  2  ♦  (94f>IUS**P/ 

«HE9E  Z  h£P8€$€NT 3  ThE  COiaoi'.iTES 
The  CIRCLE.  4 NO  ZET4  uPPRESE^TS  TiF  c-> 
090IV4TES  (*3tftFT4)  OP  T-iP  atPPOIt. 

X  a  940IUS  *  COS  ( 4NGI.I  '1.0N)  ♦  4  i  J 

Y  a  9401US  *  SIM  &NGLEvw..  . 

2  a  CMPLXCX,Y) 

2  a  2  ♦  (940!US**2)/2 
X  a  PE*L  ( Z ) 

Y  a  4  I M4G ( Z ) 

RETURN 

END 
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L' 


ti 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

to 

c 

20 

c 

30 

c 

aO 

c 

50 

C  # 
60 

C 

TO 


S'. T  I  N.f  J’  -t  (»*  .  ■■■•  ••  ) 

T *F  f  J-I C T[.r  "c  r-  ( j  ■  -  r  r  f  • 

P'jrf  t  w  e  v»L' E  'F  r-F  «  -  r) r. r [ - ■  , 

LA«)A,  GIVEN  4  V  4|_'  t  :F  ■*  \S  c -  r?  .  ’  '~‘- 
V4|-j  poir,PAMl  t,  ft  C  T '  4i  »Ty,  *  t'  i  c  <" 1  ' 

OF  LA^OAJ  t4,'04  TS  '.If  T  F  <  oi_  r; T  'I  > 

CO«Pur€c  IN  Tw?  °CUL * 4 ; < SF N  l'irFr-,3.; rp1 
4  S  <  IS.  Tw£aEF0®F .  IT  is  *.FCFS$43v  T "  ” 

OUT*  VALUER  IF  L4«C4  KNOI'O  «. 


This  SUBROUTINE  WREAKS  ThE  FuMIP-.AL  -?  ?  >_  t  T  1 .  - 

SHIP  BETWEEN  AN?  *  I“'TO  *,  E  4  =  l_  V 

LINEAR  OR  2n0  ORPER  c'.*.CTIONS  4S  *epe-o5r<tc< 

THE  V4LU?  OF  K  IS  COMPUTED  IN  ThC  a  I  \  Par. 

GRAH  4N0  IS  THEN  SE'T  TO  This  S'. '  T I  ,E .  i.  JS 
VALUE  OF  K  is  SijRJFC  TEO  to  SEVERAL  ?L  4 C f '  f  *  t 
TESTS  TO  OETERMI'‘f  T-F  approrrIatf  rf.-Iv*  r  • 

The  Lamoa-VS-K  CURVE.  L.*-w0A  ts  THE*,  fvp-/Tf' 
THROUGH  THE  APPROPRIATE  CESCRI5ING  KfL A T  I  ? •  Sh  ;  o  , 

RKl  s  -.160 
RK 2  =  -.112 
RK3  ■  n.OQ 
RK  a  a  0.06 
RK S  a  0.076 
RK6  a  0.086 
RK 7  a  O.ORUR 


IF(rk.lE.Rki) 

GO 

TO 

to 

’F(RK.LE.RK2) 

GO 

TO 

20 

IF(RK.LE.RKl) 

GO 

fO 

30 

IF(Rk.LE.pk«) 

GO 

TO 

«0 

IF(MK.LE.RKS) 

GO 

TO 

50 

IFfRK.LE.RK6) 

GO 

TO 

60 

IF(RK,GT.RK7) 

GO 

TO 

70 

RLAMOA  a  , 0  I UR» *2 

m 

(RK  -  O.o*) 

RL*VP4  a  12.- 

1  00 

,*SORT(OL JHns ) 

RE  Turn 

RLA"OA  a  C2./.0t2)  »  RK  ♦  ju.o 
RETURN 

RLAHI5A  s  (U./.ORU)  *  RK  t  2,13 

RETURN 

RLAHCU  a  (tO./,t«)  •  RK 

RETURN 


RLAHOA  a  83.33  *  RK 
RETURN 

RLAHOA  a  -  1.9  ♦  HS,  •  RK 
RETURN 

°LA*»0A  a  -6.5a  ♦  176.  *  RK 
RETURN 


RLAHR4  a  12. 

RETURN 

E*iO 
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